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1.  INTRODUCTION 


The  Subsonic  Aerodynamic  Research  Laboratory  (SARL)  wind  tunnel,  which  is  designed 
to  operate  up  to  Mach  =  0.6,  limited  to  Mach  0.5,  is  driven  by  a  20,000  HP  motor,  located 
approximately  30  ft  (9.14m)  outside  the  tunnel  exit  (Wright-Patterson  AFB,  1992).The  motor 
drives  an  approximately  2  to  1  increasing  gearbox.  The  tunnel  has  a  46  ft  by  50  ft  (14.02  m  by 
15.24  m)  inlet  area,  and  has  a  36:1  contraction  ratio  at  the  test  section.  The  tunnel  was  designed 
for  flow  visualization  and  accurate  measurements  at  low-cost  and  it  allows  for  a  360°  line-of- 
sight  of  the  model  in  the  test  section.  Turbulence  intensity  in  the  test  section  is  reported  as 
0.05%  or  below. 

Figure  1  provides  a  schematic  drawing  of  the  SARL  wind  tunnel,  and  Figure  2  provides 
a  drawing  of  the  fan  duct  and  diffuser  section  being  considered  in  this  research.  The  fan  duct  was 
designed  to  house  two-stages  of  fan  blades,  stator  vanes,  anti  swirl  vanes  and  a  nacelle.  The  fan 
was  originally  used  in  a  NASA  Langley  tunnel.  Only  the  upstream  blade  set  is  installed  in  the 
SARL.  The  nacelle  is  held  in  place  with  biconvex  struts  that  are  placed  120°  apart  from  each 
other  oriented  in  a  Y  configuration  at  the  entrance  section,  and  in  an  inverted  Y  configuration,  at 
the  aft  section  of  the  diffuser.  The  average  tunnel  velocity  at  the  entrance  to  the  fan  duct,  U,  is 
calculated  to  be  60  m/s  when  the  test  section  Mach  number  is  0.5. 

Contraction  Section 

Test  Section 


Fan  Duct 

Motor 


Figure  1  -  Schematic  drawing  with  labels  of  the  SARL  wind  tunnel  (Wright-Patterson  AFB, 

1992) 
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Figure  2  -  Drawing  of  Fan  Duct  and  Exit  Diffuser  Section  of  the  SARL  Wind  Tunnel 

(Wright-Patterson  AFB,  1992) 


Previous  analysis  of  cumulative  losses  throughout  the  tunnel  indicated  that  approximately 
30%  of  all  losses  in  the  tunnel  occurred  at  the  tunnel  exit  (Britcher,  2011).  In  that  work,  the 
geometry  of  the  entire  tunnel  was  considered  and  the  losses  were  expressed  using  a  section  total 
pressure  loss  coefficient  given  as: 


K0  = 


App  Qlocal 

Qts  Qts 


(i) 


where  AP0  is  the  change  in  total  pressure  across  a  given  location,  qts  is  the  dynamic  pressure  at 
the  test  section,  qiocai  is  the  dynamic  pressure  at  a  specified  location,  K  is  the  local  total  pressure 
loss  coefficient,  Vts  is  the  velocity  magnitude  at  the  test  section,  and  Yi0Cai  is  the  velocity 
magnitude  at  a  specified  location,  respectively.  Figure  3  provides  a  graphical  representation  of 
the  results  of  the  study,  showing  the  cumulative  total  pressure  loss  coefficient,  (£  K0 ),  versus  the 
position  within  the  tunnel. 


Figure  3  -  SARL  Cumulative  power  loss  coefficient  analysis  Britcher,  2011) 


2 

Approved  for  public  release;  distribution  unlimited. 


In  Figure  3,  the  data  points  starting  at  “Fan  Duct”  and  extending  to  “Exit”  represent  the 
portion  of  the  tunnel  that  is  being  analyzed  in  this  work.  The  results  of  the  work  by  Britcher 
(201 1)  showed  that  the  cumulative  loss  coefficient  was  about  0.245.  The  study  concluded  that 
the  circuit  losses  were  “dominated”  by  the  diffuser  and  exit  sections  of  the  tunnel.  If  the  exit  loss 
coefficient  could  be  reduced,  a  significant  reduction  in  losses  could  be  observed  for  the  entire 
tunnel,  up  to  approximately  16%  reduction  based  on  analytical  study  as  noted  by  Britcher, 

(2011). 

The  current  research  was  thus  focused  on  minimizing  the  losses  by  modifying  the 
diffuser  of  the  wind  tunnel.  In  order  to  achieve  this  goal,  several  diffuser  geometries  were 
identified  and  evaluated  for  effectiveness  in  loss  reduction.  Computational  methods  were  used 
due  to  the  very  large  size  of  the  tunnel  limiting  experimental  evaluation  of  the  models. 
Preliminary  computations  were  made  using  SolidWorks  CAD  program.  Numerical  models  of 
over  thirty  widely  varying  diffuser  geometries  were  developed  based  on  previous  knowledge  on 
diffuser  geometries  (Farokhi,  2009;  Eckert  et  al.,  1976;  Sovran  and  Klomp,  1967;  Mehta,  1977; 
Townsend,  1976;  Lefebvre,  1999;  Dixon,  2005)  and  evaluated  for  effectiveness  in  reducing 
losses  for  the  SARL  tunnel.  Literature  survey  showed  that,  previous  systematic  research 
undertaken  on  conical  diffusers  was  mainly  focused  on  studying  empty  conical  diffusers  (Fox  an 
McDonald,  1966;  Sovran  and  Klomp,  1967;  Azad,  1996;  Lefebvre,  1999;  Dixon,  2005).  These 
studies  indicate  that,  (i)  the  ratio  of  the  exit  to  entrance  areas  and  the  length  of  the  diffuser  are  the 
parameters  directly  affecting  the  pressure  recovery  performance  of  the  conical  diffusers 
(Lefebvre,  1999),  (ii)  the  pressure  recovery  becomes  independent  of  the  Reynolds  number  over 
Re  =  75000  (Fox  and  McDonald,  1966,  Sparrow  et  al.,  2009),  and  (iii)  the  diffusers  with  apex 
angle  ~7°  has  the  best  performance  (Fox  and  McDonald,  2010).  It  appears  that  research  on  the 
annular  conical  diffusers  was  focused  on  investigation  of  different  configurations  for  specific 
applications  (Johnston,  1953;  Ganesan,  1980;  Shuja  and  Habib,  1996;  Ubertini  and  Desideri, 
2000;  Cherry  et  al.,  2010;  Lo  et  al.,  2012). 

The  models  evaluated  earlier  using  the  SolidWorks  code  included  conventional  and 
unconventional  geometries,  such  as  two  diffusers,  diffusers  redirecting  the  flow  away  from  the 
ground,  diffusers  with  splitter  plates,  conical/annular  diffusers  of  different  lengths,  combination 
of  conical  diffusers,  and  constant  diameter  sections  with  varying  lengths.  Each  diffuser  model 
was  created  using  SolidWorks  CAD  program  and  CFD  solutions  were  obtained  using  the 
FloWorks  add-on  flow  analysis  tool.  The  FloWorks  solves  the  Favre  averaged  Navier-Stokes 
equations  using  the  standard  k-s  turbulence  model  (Launder  and  Spalding,  1974).  Figure  4 
illustrates  a  few  representative  examples  of  the  diffuser  geometries  considered.  Due  to  the  close 
proximity  of  the  gear  box  housing  to  the  exit  of  the  tunnel,  all  designs  had  their  lengths  restricted 
to  less  than  14  ft  (4.26  m)  since  the  diffuser  was  to  fit  between  the  tunnel  exit  and  the  gear  box 
housing  .  During  the  computations,  the  effects  of  the  blade  rotation  on  the  flow  field  were 
excluded  since  it  would  have  required  very  large  computational  resources.  The  blades  were  also 
excluded  from  the  geometry  to  study  the  effects  of  only  the  diffuser  geometry  on  the  efficiency 
calculations. 
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Figure  4  -  SARL  Example  Diffuser  Designs  Tested  using  SolidWorks 


Out  of  the  30  different  diffuser  geometries  analyzed  by  the  FloWorks  code,  4  specific 
geometries  were  selected  for  analysis  by  the  ANSYS  FLUENT  code.  These  include  the  existing 
SARL  diffuser  (baseline  model)  with  a  8°  half-apex  angle,  henceforth  termed  as  “8BT,”  and  a 
modified  version  of  it  with  a  3.5°  half-apex  angle,  henceforth  termed  as  “35BT.”  The  3.5BT  was 
modified  simply  by  adding  a  constant  area  diffuser  at  the  end,  henceforth  termed  as  “35F”  which 
was  further  modified  with  an  additional  annular  conical  diffuser  after  the  flat  section,  henceforth 
termed  as  “35FC”.  The  35FC  consisted  of  3  concentric  truncated  cones  with  included  angles  of 
60°,  34.7°,  17.2°,  and  the  entrance  area  to  the  inner  cone  was  equal  to  the  entrance  areas  between 
the  successive  cones.  The  flat  and  the  conical  sections  had  the  same  length.  The  35F  and  the 
35FC  configurations  were  selected  because  they  produced  greatest  reduction  in  pressure  losses 
and  about  6%  improvement  in  the  entire  tunnel  efficiency  when  analyzed  with  the  FloWorks 
code.  These  4  diffuser  configurations  are  shown  in  Figure  5. 
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Figure  5  -  Sectional  views  of  the  4  diffuser  configurations;  a)  8BT,  b)  35BT,  c)  35F,  d)  35FC 

For  each  of  these  models,  both  uniform  and  fully  developed  inlet  flow  conditions  were 
used.  Computations  were  made  on  full  scale  models  of  the  existing  SARL  tunnel.  Thus 
computations  were  performed  for  a  total  of  8  different  configurations.  Table  1  presents  the 
summary  matrix  of  the  various  computations  performed  in  this  study.  The  Reynolds  number  for 
these  computations  is  based  on  the  average  inlet  velocity  and  the  duct  diameter  at  the  inlet 
section  for  air  at  15°C.  The  Reynolds  number  for  the  full  scale  models  were  taken  as  1.647><107. 

Table  1  -  Full  Matrix  of  Computations  by  Model  and  Boundary  Conditions 


3D  Full  Size 
(Re  =  1.647xl07) 

Grid  independence  study 

8BT,  35FC 

Uniform  inlet  flow 

8BT,  35BT,  35F,  35FC 

Fully  developed  inlet  flow 

8BT,  35BT,  35F,  35FC 
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2.  NUMERICAL  PROCEDURE  AND  PROBLEM  SETUP 

The  Reynolds-Averaged  Navier-Stokes  (RANS)  equations  were  solved  using  the  ANSYS 
FLUENT  commercial  code  with  the  standard  k-s  turbulence  model  (Launder  and  Spalding, 

1974).  The  governing  equations  and  turbulence  model  constants  are  available  in  FLUENT  user 
manual  and  are  not  repeated  here.  The  flow  domain  is  divided  into  many  small  finite  volumes 
with  an  unstructured  mesh.  A  collocated  arrangement  for  the  placement  of  the  flow  variables  is 
used  in  the  mesh  system.  The  conservation  equations  are  integrated  over  each  of  the  finite 
volume  to  yield  sets  of  linear  algebraic  equations.  These  sets  of  linear  algebraic  equations  are 
then  solved  sequentially  using  an  iterative  method.  The  SIMPLE  (Patankar,  1980)  algorithm  is 
used  for  the  pressure -velocity  coupling.  The  convective  fluxes  have  been  calculated  using  the 
second  order  upwind  scheme  while  the  diffusive  fluxes  have  been  calculated  using  the  central 
difference  scheme.  The  convergence  is  assumed  when  the  value  of  the  scaled  residual  of 
continuity  and  momentum  equations  is  less  than  10~6.  Convergence  has  also  been  monitored  by 
plotting  the  drag  on  the  nacelle  surface  until  the  variation  of  the  drag  leveled  off  with  iteration. 

2.1  Geometry  and  Mesh  Design. 

The  flow  configuration  geometry  was  created  using  the  SolidWorks  code  which  was 
imported  into  ANSYS  FLUENT  for  mesh  generation.  The  computational  mesh  was  generated 
using  the  meshing  code  integrated  in  the  ANSYS  package.  The  solid  bodies,  such  as  the  outer 
diffuser,  fan  duct  wall,  engine  nacelle,  and  support  spars  had  a  “face  sizing”  mesh  control 
applied  to  these  surfaces,  with  a  slow  smoothing  method  applied  to  the  domain  from  these 
locations.  This  created  an  unstructured  grid  that  had  the  highest  mesh  refinement  near  the  actual 
diffuser  geometry,  which  slowly  became  coarser  farther  away  from  these  surfaces.  The  flow 
domain  consists  of  two  parts;  the  base  tunnel  fan-duct  and  diffuser  from  the  inlet  to  the  exit  (the 
region  between  the  outer  casing  and  nacelle  surface,  the  downstream  diffuser  section,  and  the 
upstream  inlet  duct  section)  and,  the  downstream  open  rectangular  wake  region  as  shown  in 
Figure  6.  A  close  up  view  of  the  sample  base  configuration  surface  mesh  is  shown  in  Figure  7. 
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Figure  6  -  Sample  Mesh  of  Base  Diffuser  Geometry  for  Entire  Computational  Domain  for  8 

Base  Configuration  in  ANSYS  FLUENT. 


Y 


Figure  7  -  Sample  Mesh  of  Base  Diffuser  Geometry  showing  Surface  Meshing  for  8  Base 

Configuration  in  ANSYS  FLUENT. 

The  three-dimensional  domain  required  a  large  number  of  cells  (of  the  order  of  a  few  million)  to 
fill. 
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2.2  Boundary  Conditions 

For  all  computations  the  no-slip  condition  was  imposed  on  all  solid  walls  or  surfaces 
including  the  nacelle  surface,  tunnel  casing  wall,  on  the  strut  rod  surfaces,  etc.  Constant  pressure 
outlet  condition  was  imposed  on  all  open  boundaries  of  the  downstream  rectangular  wake  region 
where  the  values  of  the  velocity  components  are  extrapolated  from  inside.  As  mentioned  earlier, 
two  types  of  inlet  conditions  were  used  in  the  computations;  (/')  uniform  axial  velocity 
compatible  to  the  corresponding  Reynolds  number  with  other  velocity  components  being  zero 
and  (ii)  a  fully  developed  mean  axial  velocity  component  with  average  velocity  compatible  to  the 
Reynolds  number.  The  uniform  inlet  velocity  condition  is  straight  forward  to  implement.  For  the 
fully  developed  inlet  velocity  condition,  the  profile  was  obtained  through  separate  turbulent  flow 
simulation  in  a  circular  pipe  of  length  greater  than  the  entry  length,  which  is  given  by 
4.4(D)(Re1/6 ),  where  D  is  the  pipe  diameter  and  Re  is  the  Reynolds  number  (given  in  Table  1) 
based  on  the  average  velocity.  A  uniform  velocity  of  60m/s  was  input  as  the  velocity  entering  the 
pipe  and  the  profile  at  the  exit  of  this  pipe  was  specified  as  the  tunnel  section  inlet  velocity 
condition  through  the  user  defined  function  feature  of  the  FLUENT  code. 

2.3  Grid  Independence  Study 

In  order  to  obtain  mesh  independent  results,  systematic  mesh  independence  study  was 
conducted.  Starting  with  a  basic  coarse  mesh  with  few  hundred  thousand  cells,  the  mesh  was 
successively  refined  and  the  computed  shear  stress  distributions  along  the  nacelle  top  surface  as 
well  as  the  total  pressure  profiles  at  the  diffuser  exit  along  the  y-mid-pane  were  compared.  The 
mesh  was  considered  optimally  refined  when  the  plots  of  the  shear  stress  and  total  pressure 
distribution  did  not  change  appreciably  between  successive  mesh  refinements. 

Figures  8  and  9  provide  examples  of  the  results  of  grid  independence  studies  for  the 
corresponding  geometry  mentioned  in  the  plot  title.  Once  grid  independence  was  established  for 
the  original  diffuser  geometry,  the  same  mesh  settings  were  applied  to  the  other  geometries,  as 
the  changes  were  assumed  relatively  minor.  One  exception  to  this  was  the  diffuser  geometry  that 
made  use  of  conical  sections  to  split  the  flow  at  the  exit  of  the  diffuser  (35FC  case).  A  separate 
grid  independence  study  was  made  for  this  geometry. 
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Figure  8  -  Variation  of  total  pressure  at  the  diffuser  section  exit  plane,  along  the  central 
line  in  the  y-direction,  with  mesh  refinement.  Numbers  at  the  bottom  represent  total 

number  of  cells  for  each  mesh. 
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Figure  9  -  Variation  of  shear  stress  along  engine  nacelle  top  surface  with  mesh  refinement. 
Numbers  at  the  bottom  represent  total  number  of  cells  for  each  mesh 
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3.  RESULTS  AND  DISCUSSION 

For  the  SARL  tunnel,  the  head  loss  is  considered  in  two  parts:  the  head  loss  within  the 
fan  duct  and  diffuser,  and  the  head  loss  of  the  flow  exiting  the  diffuser.  Figure  10  represents  the 
locations  for  which  these  losses  are  associated.  Here  location  1  and  2  represent  the  entrance  and 
exit  sections  of  the  diffuser  while  location  3  represents  a  section  outside  the  exit  of  the  diffuser 
which  is  sufficiently  far  away  so  that  the  velocity  at  this  location  is  essentially  zero.  The  total 
head  loss  for  sections  1-2  and  2-3,  can  be  expressed  with  the  following  equations: 


,  _  (Pi-P2)+\p(.v?-v£) 

hLi_2  — 


pa 


hL 


2-3 


(2) 

(3) 


h 


Ltotal 


=  hLl_  2  +  hLz_  3 


(4) 


Figure  10  -  General  locations  for  power  loss  calculations  relative  to  the  SARL  wind  tunnel 

fan  duct  and  diffuser  sections. 

Since  the  head  losses  are  directly  related  to  the  change  in  the  total  pressure  from  the  inlet 
to  the  exit  of  the  diffuser,  the  total  pressure  distributions  both  in  the  fan  duct  and  the  diffuser 
sections  are  investigated.  The  total  pressure  distribution  for  the  original  configuration,  8BT,  with 
uniform  inlet  velocity  condition  is  shown  in  Figure  11.  The  flooded  contours  of  the  total 
pressure  distribution,  at  the  diffuser  exit  section  and  on  a  longitudinal  plane  in  the  axial  direction, 
are  shown  on  the  left  and  right,  respectively,  in  this  figure.  The  longitudinal  flow  direction  in  this 
and  other  subsequent  figures  is  from  left  to  right.  Additional  contours  of  velocity  magnitude 
distribution  through  the  geometry  for  the  same  case  are  provided  in  Figure  12  to  visualize 
potential  locations  of  flow  separation  and  diffusion. 
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Figure  11  -  Total  pressure  distribution  through  8BT  -  uniform  inlet  flow. 


Figure  12  -Velocity  magnitude  distribution  through  8BT  -  uniform  inlet  flow. 

It  can  be  seen  that  the  flow  does  begin  to  slow  significantly  along  the  8°  angle  diffuser 
section,  and  the  velocity  distribution  at  the  exit  is  generally  uniform  and  axisymmetric.  To 
compare  this  to  the  fully  developed  flow  condition  at  the  inlet,  the  distribution  of  the  velocity 
through  the  model  has  been  provided  in  Figure  13. 
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Figure  13  -  Velocity  magnitude  distribution  through  8BT  -  fully  developed  inlet  flow. 

While  the  velocity  distribution  within  the  fan  duct  and  diffuser  does  differ  slightly  from 
that  of  the  uniform  flow  case,  the  velocity  distribution  at  the  exit  does  have  one  region  of  notable 
difference  from  the  uniform  inlet  flow  case.  The  velocity  distribution  at  the  exit  of  the  8BT  for  a 
fully  developed  inlet  flow  is  generally  uniform,  but  shows  a  region  of  increased  velocity  along 
the  bottom  of  the  exit  between  120°  and  240°  counterclockwise  from  the  top  of  the  exit  section. 
This  altered  velocity  distribution  could  be  attributed  to  the  flow  remaining  attached  a  greater 
distance  for  fully  developed  flow  than  for  uniform  flow  within  the  diffuser,  and  to  the 
interference  of  a  support  spar  directly  immediately  upstream  of  the  exit. 

Next,  the  35BT  geometry  is  considered.  Figure  14  represents  the  total  pressure 
distribution  through  the  model  for  the  uniform  inlet  flow  case,  while  the  Figure  15  represents  the 
velocity  magnitude  distribution  through  the  geometry. 


Figure  14  -  Total  pressure  distribution  through  35BT  -  uniform  inlet  flow. 
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Figure  15  -  Velocity  magnitude  distribution  through  35BT  -  uniform  inlet  flow. 

Compared  to  the  distributions  of  total  pressure  and  velocity  magnitude  observed  in  the 
8BT  model,  it  seems  that  there  is  a  greater  degree  of  diffusion  occurring  in  the  flow  near  the  exit 
of  the  35BT  diffuser  that  is  likely  to  lead  to  greater  losses.  To  visualize  and  compare  the 
differences  between  the  uniform  inlet  flow  and  fully  developed  flow  cases,  the  velocity 
distribution  through  the  full  size  35BT  geometry  for  fully  developed  inlet  flow  is  shown  in 
Figure  16. 


Figure  16  -  Velocity  magnitude  distribution  through  35BT  -  fully  developed  inlet  flow. 

For  both  the  uniform  and  fully  developed  flow  cases,  the  flow  travelling  through  the 
35BT  geometry  encounters  the  same  diffusion  as  it  nears  the  exit  of  the  diffuser. 

A  similar  comparison  of  the  total  pressure  and  velocity  magnitude  distributions  of  the 
35F,  and  35FC  models  is  presented  for  comparison  to  the  baseline  geometry.  Figure  17 
represents  the  total  pressure  distributions,  while  Figure  18  represents  the  velocity  magnitude 
distributions  through  the  full  size  35F  diffuser  model  under  uniform  inlet  flow  conditions. 
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Figure  18  -  Velocity  magnitude  distribution  through  35F  -  uniform  inlet  flow. 

The  figures  indicate  that  the  velocity  of  the  flow  exiting  the  35F  full  size  model  is  lower 
than  that  of  the  baseline  geometry  (35BT).  In  addition,  the  flow  can  be  observed  to  slow 
following  the  surface  of  the  engine  nacelle  near  the  exit  of  the  diffuser,  compared  to  that  of  the 
35BT. 

The  Figure  19  shows  the  velocity  distribution  through  the  35F  under  fully-developed 
inlet  conditions.  Velocity  distribution  is  similar  to  the  distribution  obtained  with  the  uniform  inlet 
flow  case.  The  flow  velocity  peaks  in  a  doughnut  section  closer  to  the  mid  section  of  the  duct, 
the  flow  near  the  nacelle  has  decreased  velocity  values. 
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Figure  19  -  Velocity  magnitude  distribution  through  35F  -  fully  developed  inlet  flow 


Figures  20  and  21  represent  the  total  pressure  and  velocity  distributions,  respectively,  of 
the  full  size  35FC  diffuser  geometry. 


Figure  20  -  Total  pressure  distribution  through  35FC  -  uniform  inlet  flow. 
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Figure  21  -  Velocity  magnitude  distribution  through  35FC  -  uniform  inlet  flow. 

Figures  indicate  that  the  total  pressure  exiting  the  diffuser  is  essentially  at  atmospheric 
conditions  and  that  the  flow  through  the  fan  duct  section  is  free  of  separation.  However,  within 
the  conical  region  of  the  diffuser,  flow  separation  occurs  along  the  surface  of  the  cone  splitting 
the  flow.  The  conical  sections  that  are  radially  further  from  the  center  of  the  diffuser  (between 
the  outer  and  middle  cones)  have  more  separation  then  the  middle  sections  (inner  cone  and 
between  the  inner  and  middle  cones).  In  addition,  the  flow  that  is  radially  closest  to  the  center  of 
the  diffuser  is  turned  by  the  inner  conical  section,  increasing  the  spreading  rate  of  the  jet  exiting 
the  tunnel.  Considering  the  total  pressure  distribution  at  the  exit  of  this  diffuser  geometry,  it  is 
expected  that  the  35FC  model  will  offer  the  greatest  improvement  in  the  reduction  of  losses. 

The  flow  velocity  magnitude  obtained  with  fully  developed  inlet  flow  conditions  is 
shown  in  Figure  22.  Figure  shows  that  the  velocity  distribution  is  very  similar  to  the  one 
obtained  with  uniform  inlet  flow  conditions  but  the  velocity  magnitude  is  observed  to  be  smaller 
throughout  the  fan-duct,  diffuser  and  at  the  exit. 


Figure  22  -  Velocity  magnitude  distribution  through  35FC  -  fully  developed  inlet  flow. 

16 

Approved  for  public  release;  distribution  unlimited. 


The  total  head  losses  associated  with  the  various  diffuser  configurations  are  computed 
using  Equations  (2)  -  (4)  and  are  presented  in  Table  2  for  each  model  considered  in  this  study. 

Table  2  -  Calculated  head  losses  (m)  at  standard  pressure  and  temperature  for  various 


difl 

fuser  geometries. 

Diffuser 

Uniform 

Fully  Developed 

Configuration 

Inlet  Flow 

Inlet  Flow 

8BT 

113.04 

106.30 

35BT 

153.80 

136.12 

35F 

108.30 

94.54 

35FC 

96.24 

84.12 

Results  indicate  that  a  fully  developed  inlet  flow  leads  to  a  reduction  for  all  the  diffuser  models, 
and  that  the  35FC  is  the  best  working  model. 

A  comparison  of  the  percentage  improvement  of  each  model  relative  to  the  base  line 
(8BT)  configuration  is  compiled  into  Table  3. 

Table  3  -  Percentage  improvement  in  head  losses  compared  to  the  base  line  (8BT) 
_ configuration. _ 


Diffuser 

Configuration 

Uniform 

Inlet  Flow 

Fully  Developed 
Inlet  Flow 

35BT 

-36.05 

-28.05 

35F 

4.20 

11.07 

35FC 

14.87 

20.86 

It  can  be  seen  that  for  the  full  size  models,  the  35FC  diffuser  geometry  offers  the  greatest 
improvement  to  reducing  losses  with  a  range  between  14.87%  for  the  uniform  inlet  flow  case  to 
20.86%  for  the  fully  developed  flow  case.  Considering  the  entire  tunnel,  this  translates  to 
approximately  a  5.2%  to  7.3%  improvement.  The  35F  geometry  also  offers  an  improvement 
range  of  4.20%  to  1 1.07%,  which  translates  to  1.47%  to  3.87%  range  for  the  entire  tunnel.  Thus, 
from  the  numerical  analysis  presented  here  it  is  concluded  that  the  35FC  diffuser  geometry  yields 
maximum  pressure  recovery  in  the  full  size  tunnel. 

3.1  Further  Discussion 

As  discussed  by  Lefebvre  (1999)  based  on  the  research  conducted  by  Sovran  and  Klomp 
(1967)  the  aspect  ratio,  AR  =  A2/A1,  and  the  dimensionless  diffuser  length,  DL  =  L/R1  of  the 
conical  or  annular  conical  diffusers  define  the  pressure  recovery  of  the  diffuser  rather  than  the 
geometry  of  the  annular  sections.  In  these  equations,  At,  and  A2  are  the  entrance  and  the  exit 
areas  of  the  diffuser,  L  is  the  length  of  the  diffuser  and,  the  R1  =  -J  A1/n  is  the  radius  of  the 
diffuser  at  the  entrance.  The  performance  curves  generated  by  Sovran  and  Klomp  (1967)  are 
shown  in  Figure  23.  The  parameters  pertaining  to  the  current  diffuser  geometries  were 
calculated  and  are  presented  in  Table  4.  Parameters  were  calculated  using  two  separate 
arguments,  in  the  first  case  calculations  were  done  as  if  the  diffuser  section  could  be  assumed  as 
an  empty  conical  section  and  thus  the  nacelle  was  neglected,  in  the  second  case  cross  sections 
that  were  blocked  by  the  nacelle  were  omitted.  Also,  in  the  second  case  the  entrance  radius  was 
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calculated  as  the  radius  of  a  circle  that  has  an  area  equal  to  that  of  the  unblocked  area.  The 
coordinates  of  the  geometries  are  shown  in  Figure  23. 


Dimensionless  diffuser  length,  DL  =  L/R± 


Figure  23  -  Performance  chart  for  conical  diffusers,  (adapted  from  Sovran  and  Klomp, 

1967). 

The  Cp*  line  in  Figure  23  is  the  locus  of  points  which  define  the  diffuser  area  ratio 
producing  maximum  pressure  recovery  in  a  prescribed  non-dimensional  diffuser  length.  The 
Cp**  line  is  the  locus  of  points  which  define  the  diffuser  non-dimensional  length  producing 
maximum  pressure  recovery  at  a  prescribed  area  ratio.  Cp  is  the  pressure  recovery  coefficient 

based  on  the  bulk  inlet  velocity.  The  symbols  ""I  "and  denotes  the  empty  conical  and  annular 

conical  cases,  respectively. 
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Table  4-  Calculated  head  losses  (m)  at  standard  pressure  and  temperature  for  various 

diffuser  geometries. 


geometry 

empty 

conical 

annular 

conical 

AR 

DL 

AR 

DL 

8BT 

1.34 

1.13 

1.83 

1.43 

35BT 

1.15 

1.13 

1.51 

1.43 

35F 

1.15 

2.81 

1.83 

3.56 

35FC 

2.39 

2.81 

3.54 

3.57 

The  parameters  were  computed  in  two  different  ways  since  it  was  observed  that  the  8BT 
parameters  indicated  to  a  point  very  close  to  the  best  pressure  recovery  line  once  empty  conical 
argument  is  made,  however  the  parameters  indicate  to  a  point  above  the  0=10°  curve  for  the 
annular  conical  argument  indicating  to  a  possibility  that  a  different  diffuser  could  be  used  to 
improve  the  performance  of  the  SARL.  The  correct  use  of  the  chart  requires  using  only  the  areas 
that  allow  fluid  flow.  The  35BT  parameters  indicate  that  the  annular  conical  calculations  result  in 
a  point  near  the  best  pressure  recovery  line,  although  CFD  results  showed  that  this  geometry 
performed  poorly.  The  parameters  calculated  for  the  35F  differ  drastically  from  each  other.  The 
annular  conical  argument  indicates  to  the  fact  that  this  diffuser  is  the  most  likely  candidate  since 
the  geometry  of  the  diffuser  is  easy  to  manufacture  and  it  is  the  second  best  performing 
geometry.  The  35FC  diffuser  geometry  that  resulted  in  best  performance  in  CFD  calculations 
resides  away  from  the  best  performance  curve  shown  by  Sovran  and  Klomp  (1967).  A  similar 
geometry  was  discussed  by  Cheng  (1992)  as  a  means  to  increase  the  pressure  recovery  while 
reducing  the  diffuser  length.  Cheng  had  argued  that  different  conical  sections  should  have  a 
common  apex  point,  and  that  the  flow  can  be  considered  as  spherically  diffusing.  The  35FC 
geometry  differs  than  the  one  discussed  by  Cheng  (1967)  since  the  conical  annular  sections  are 
same  length,  ft  is  believed  that  the  conical  sections  distribute  and  force  the  flow  velocity  to  slow 
down  although  they  also  cause  separation  in  the  conical  sections.  For  the  application  studied  in 
the  current  paper,  the  separation  closer  to  the  exit  is  not  important  since  there  are  no  tunnel 
sections  beyond  this  point.  However,  separation  may  cause  increased  noise  levels,  which  will  be 
studied  in  a  future  study. 
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4.  CONCLUSIONS 

The  primary  goal  of  this  research  was  to  identify  the  diffuser  geometry  that  would  results 
in  greatest  reduction  in  pressure  losses  for  the  Subsonic  Aerodynamic  Research  Laboratory 
(SARL)  wind  tunnel  at  the  Air  Force  Research  Laboratory.  Computational  results  of  two  diffuser 
geometries  were  compared  to  each  other  and  to  the  existing  diffuser  in  providing  the  maximum 
pressure  recovery.  Flow  field  solutions  were  obtained  at  60  m/s  average  inlet  velocity  for  both 
uniform  inlet  turbulent  flow  and  fully  developed  inlet  turbulent  flow. 

Results  indicate  that  the  35FC  diffuser  geometry  offered  the  greatest  reduction  in 
pressure  losses,  yielding  14.87%  to  20.86%  improvement  over  the  existing  SARL  tunnel  diffuser 
for  the  full  scale  model.  Considering  the  improvement  for  the  entire  SARL  tunnel,  this  range  is 
5.20%  to  7.30%.  The  35F  geometry  also  offered  improvement  in  the  reduction  of  losses, 
resulting  in  4.20%  to  1 1 .07%  improvement  over  the  existing  SARL  diffuser.  For  the  entire 
SARL  tunnel,  this  range  is  approximately  1.47%  to  3.77%  improvement  in  pressure  recovery. 
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